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\D . ABSTRACT 

' We report on an ambitious multi-site campaign aimed at detecting stellar variability, 

fH \ particularly solar-like oscillations, in the red giant stars in the open cluster M67 (NGC 

2682). During the six- week observing run, which comprised 164 telescope nights, we 
used nine 0.6-m to 2.1-m class telescopes located around the world to obtain uninter- 
rupted time-series photometry. We outline here the data acquisition and reduction, 
with emphasis on the optimisation of the signal-to-noise of the low amplitude (50- 
500 /imag) solar-like oscillations. This includes a new and efficient method for obtain- 
ing the linearity profile of the CCD response at ultra high precision 10 parts per 
million). The noise in the final time series is 0.50 mmag per minute integration for the 
best site, while the noise in the Fourier spectrum of all sites combined is 20 /zmag. In 
addition to the red giant stars, this data set proves to be very valuable for studying 
' high-amplitude variable stars such as eclipsing binaries, W UMa systems and 5 Scuti 

stars. 

Key words: Stars: red giants - Stars: oscillations - Stars: variables: delta Scuti - 
Stars: variables: W UMa - Stars: blue stragglers - Open clusters: individual (NGC 
2682, M67) - Techniques: photometric. 
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1 INTRODUCTION 

Asteroseismology of stellar clusters is potentially a power- 
ful tool. The assumption of a common age, distance, and 
chemical composition provides stringent constraints on each 
cluster member, which significantly improves the astero- 
seismic output fCough fc Novotnvl Il99jh . Hence, detect- 
ing oscillations in cluster stars in a range of evolutionary 
states holds promise of providing new tests of stellar evolu- 
tion theory. Driven by this great potential, several studies 
have been aimed at detecting solar-like oscilla tions in the 
open cluster M67 f Gilliland & Brown 1988; Gillil and et"al] 
Il99ll: iGilliland fcljrowni il992a : iGilliland et al.i il993f) and 
in the globular cluster M4 (Frandsen et al. in prep.). The 
most ambitious campaign was reported by IGilliland et alJ 



l|l993r . who used seven 2.5-m to 5-m class telescopes dur- 
ing one week in a global photometric network to target 
11 turn-off stars in M67. Despite these efforts, they did 
not claim unambiguous detection of oscillations. However, 
one of their conclusions was that oscillations should be de- 
tectable in the more evolved red giant stars due to higher 
expected oscillation amplitudes. However, the oscillation pe- 
riods of up to several hours and expected frequency sep- 
arations of a few 10~^ Hertz would require a time base 
of roughly one month on these stars. Recent month-long 
studies using single- or dual-site high-precision radial ve- 
locity measurements (cr ~ 2m/s) on bright field stars have 
clearly demonstrat ed that solar-like oscillati ons are present 
in red giant star s ( Fran dsen et aljEoO l lBarban et al.ll20o3 : 
Ide Ridder et alj|2006i'l . However, due to non-continuous cov- 



2 D. Stello et al. 




(B-V) 

Figure 1. Col our-Magnitude diagram of the open cluster M67 
(photometry bv lMontgomerv et alj ll993l) '). The red giant target 
stars are indicated with filled symbols. The numbers correspond 
to those indicated in Fig. The solid line is an isochrone ((m — 
M) = 9.7mag, Age = 4.0 Gyr. Z = 0.0198 an d Y = 0.2734) from 
the BaSTI database jPietrinferni et al.ll2004h . 
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Figure 2. Time series of star No. 10 for all sites (after remov- 
ing outliers and correcting for colour-dependent extinction; see 
Sect. 1^^ . The site abbreviations are explained in Table HI 

Fourier spectra of these stars will be presented by Stello et 
al. (in prep.). 



erage these data suffered badly from aliasing in the Fourier 
spectrum, which complicated the detailed frequency analy- 
sis. Multi-sit e or space observations are therefore required 
JStello et alj |2006). Such observations will hopefully soon 
become available for red giant stars in the field from the 
MOST, COROT and Kepler missions. However, after the 
cancellation of the ESA Eddington mission, no current 
or planned space project will measure stellar oscillations 
in cluster stars. Using high-precision spectrographs from 
ground to measure radial velocities in red giants is not pos- 
sible due to the lack of a global network that can achieve 
high-precision velocities on an ensemble of relatively faint 
cluster stars. Hence, the only feasible approach is ground- 
based photometry. 

In this investigation we again target M67 using multi- 
site photometric observations. Unlike the previous stud- 
ies on this cluster, our primary targets are the red giant 
stars (see Fig. [Jl. Ex trapolating the L/Af -scaling relation 
ijKieldsen fc Beddind [l995') predicts the amplitude of these 
stars to be in the range 50-500/imag. Although very low, 
these amplitudes are sig nificantly higher than for the turn- 
off stars targeted by, e.g.. lGilliland et alJ (Il99.?l) . In addition 
to the red giants, more than 300 cluster stars were observed 
during the campaign. Many are high-amplitude variables 
such as W UMa systems, 5 Scuti stars and eclipsing binaries, 
some of which are already known. This campaign provides 
a unique data set to investigate these stars as well (Bruntt 
et al., in prep.). 

With emphasis on the low-amplitude red giant stars, 
the main purpose of this paper is to describe the optimisa- 
tion of time-series data towards achieving the highest pos- 
sible signal-to-noise in the Fourier spectrum (in amplitude). 
Further discussion on the extraction of p-modes from the 



2 OBSERVATIONS 

We observed the open cluster M67 from 6 January to 17 
February 2004 using nine telescopes (0.6-m to 2.1-m class) 
in a global multi-site network. The sites were distributed in 
longitude to allow continuous time series photometry dur- 
ing the six-week observing program. We were allocated 164 
nights of telescope time which, due to bad weather, yielded 
about 100 clear nights (see Fig. |5|and Table 0. In the first 
18 days we observed 34% of the time and in the following 
three weeks the coverage was 80%. For the entire campaign 
(43 days) the coverage was 56%. 

The participating telescopes and detectors had different 
properties and the data sets are therefore rather diverse in 
terms of field-of-view (FOV), cadence and noise properties. 
A summary of the observations and the instrument charac- 
teristics for each site is given in Table We indicate the 
smallest and largest FOV in Fig. |3 The red giant stars are 
indicated as well. 

Observations at each telescope were planned to opti- 
mise the signal-to-noise for solar-like oscillations in the red 
giant stars. We did that by calculating both the noise and 
the expected oscillation amplitudes (the signal) in differ- 
ent photometric filters. The amplitudes were estimated from 
iKieldsen fc Beddind (Il995l Eq. 5) and the noise was es- 
timated by photon counting statistics. These calculations 
showed that the Johnson B and V filters were favourable, 
but an on-site test was required to establish which of these 
was superior at each telescope. The observers therefore chose 
filters based on an initial test at the beginning of the first 
night. All sites except Kitt Peak chose the V filter. No phase 
change is seen in the solar oscillations between observa- 
tions obtained in different filters over the range 400-700 nm 
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J Jimenez et alJll999l) . We therefore expect the same adia- 
batic behaviour for high-order solar-like oscillations in other 
stars as well. Hence, combining data based on different fil- 
ters can be done after a simple rescaling of the amplitude, 
and corresponding adjustment of the weights to preserve the 
signal-to-noise. 

To obtain a noise level which was essentially limited 
by scintillation and photon noise, it was important to avoid 
drift on the CCD of the stellar field. The aim was to have 
each star confined within a few pixels. Not all sites had au- 
toguiding systems and as a result we found very different 
drift characteristics from site to site (see Fig. |3J . The im- 
ages were defocused to obtain a higher duty cycle but we 
avoided crowding. 

The exposure time at each site was adjusted to have 
star No. 4 just below the saturation limit, which provided a 
safety margin for the large group of clump stars that were 
0.3 mag fainter (see Fig.0. However, the two brighter stars 
(No. 3 and 11) were therefore often saturated. Due to their 
expected long oscillation periods (Stello et al., in prep.), on 
time scales similar to typical instrumental drift, the results 
on these stars were likely to provide only limited scientific 
output. Keeping these two stars above the saturation limit 
resulted in lower noise for the stars at the base of the red 
giant branch (stars 6, 13, 12, 14) which were more likely to 
produce useful results. 



3 CALIBRATION 

We calibrated each CCD image using four steps: 

(i) overscan subtraction (not all CCDs had an overscan 
region), 

(ii) subtraction of bias (the bias levels were stable enough 
to use a single master bias image for each CCD), 

(iii) correction for non-linearity (see Sect. 13.11 . 

(iv) fiat-fielding to correct for pixel-to-pixel variations in 
the quantum efficiency (we used one master fiat field for 
each CCD; for Kitt Peak and the RCC this was based on 
dome fiats, while we used sky fiats for all other sites). 

We found that the dark current was negligible compared to 
the read-out-noise for all sites and it was therefore ignored. 
These four steps were standard except the non-linearity cor- 
rection, which is described in more detail in the following 
section. 



3.1 CCD linearity calibrations 

In this project a few target stars were relatively close to 
the CCD saturation limit, at fiux levels for which the non- 
linear response of the CCD gain could be significant. It is 
important to correct for these gain variations to attain the 
high photometric precision required by this project. 

Linearity at high fiux levels was investigated for all 
CCDs using a "classical" linearity test. We used an ap - 
proach similar to that described bv lGilliland et aTl il993l) . 
The method measures relative variations in the CCD am- 
pHfier gain, rather than absolute calibration in terms of 
e~/ADU. Flat-field images were taken sequentially with in- 
creasing exposure time, interleaved with reference images; 
e.g. 3s, 10s, 3s, 10s,- ••, 3s, then followed by 3s, 20s, 3s, 
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Figure 4. The position on the CCD of star No. 10 relative to 
a reference frame. The insets show the inner 20 by 20 pixels. 
Although autoguiding is good at LaS, instrument rotation intro- 
duced drift during the observing run for stars far from the rotation 
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Table 1. Summary of observations. 



Site*^ 


Telescope 


Filter 


FOV 


Image 


Observing 


Nexp 


Median 


Exp. 


Duty 


#nights 


#nights 




tx[JK:L L Ul C 






scale 
















SSOi 


1.0m 


V 


14^0 


0^'38/pix 


43.0h 


2205 


62s 


30s 


48% 


16 


9 


SSO2 


1.0m 


V 


14.'0 


O'.'60/pix 


51. Oh 


1157 


144s 


30s 


21% 


17 


10 


SOAO 


0.6m 


V 


20.'5 


0y60/pix 


33.6h 


467 


240s 


120s 


50% 


17 


7 


SAAO 


1.0m 


V 


6^0 


0"31/pix 


112.9h 


2595 


149s 


80s 


54% 


28 


22 


RCC 


1.0m 


V 


T.O 


0('29/pix 


23.6h 


722 


89s 


70s 


79% 


14 


5 


Sch 


0.6m 


V 


17.'0 


l^lO/pix 


31. 3h 


1584 


53s 


35s 


66% 


16 


6 


LaS 


1.5m 


V 


13^5 


O'.'39/pix 


lOQ.Oh 


3945 


90s 


24s 


27% 


22 


18 


LOAO 


1.0m 


V 


22^5 


O'.'66/pix 


41.8h 


2886 


46s 


12s 


26% 


14 


7 


Kitt 


2.1m 


B 


lO'.O 


0'.'30/pix 


75.5h 


1563 


172s 


52s 


30% 


11 


9 


Lag 


1.0m 


V 


14^0 


O'.'41/pix 


46.4h 


1320 


114s 


20s 


18% 


9 


7 


Total 










561. Ih 


18444 








164 


100 



" Site abbreviations and observers are: SSOi (Wide Field Imager at Siding Spring Observatory, Australia, Z.E.D., 
D.S., A.P.J, and L.L.K.); SSO2 (Imager at Siding Spring, A.P.J. , J.N. and D.S.); SOAO (Sobaeksan Optical Astronomy 
Observatory, Korea, S.-L.K., J.-A.L. and C.-U.L.); SAAO (South Africa Astronomical Observatory, T.A. and H.R.J. ); 
RCC (Ritchey-Chretien-Coude at Piszkesteto, Konkoly Observatory, Hungary, J.N.); Sch (Schmidt at Piszkesteto, Z.C. 
and R.S.); LaS (La Silla Observatory, Chile, H.B.); LOAO (Mt. Lemmon Optical Astronomy Observatory, Arizona, Y.B.K. 
and J.-R.K.); Kitt (Kitt Peak National Observatory, Arizona, R.L.G.); Lag (Mt. Laguna Observatory, California, C.S. 
and M.Y.B.). 
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Figure 3. M67 ficld-of-vicw (FOV). SAAO with the smallest FOV and LOAO with the largest are indicated with dashed squares. The 
red giant stars are marked with circles. Source: STScI Digitized Sky Survey. 
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Figure 6. Linearity flat field of tlie longest exposure, 410 s (see 
text). Tfie fiat plateau at high row numbers is due to the digital 
saturation of the CCD. 



Figure 5. Classical linearity tests for all CCDs. For SSOi, SSO2 
and SAAO two separate tests have been merged. Error bars are 
plotted as 3xrms to make them visible. The solid lines are poly- 
nomial fits used to correct the data for non-linearity. 



20s,- • •, 3s, until the final series, in wiiicii the longer expo- 
sures were almost saturated. Tiie reference images allowed 
instabilities of the light source to be measured and removed. 
In some cases, however, the flat-field lamp varied on time 
scales too short to be sampled and a correction could there- 
fore not be made. The mean counts in the flat-field image, 
scaled according to the exposure time, were plotted versus 
the mean counts. Images with the same exposure time were 
grouped to form a single point, with an uncertainty equal to 
the group rms. In Fig.|S]we show the results of the linearity 
tests for all CCDs. Saturation occurred at 65536 ADU except 
for the Schmidt, where it was at 16384 ADU. From Fig. H 
(upper right panel) we see that non-linearity can introduce 
variability in the stellar time series up to several percent 
in non-photometric conditions (variable atmospheric trans- 
parency). Using ensemble photometry (Sect. 0J will, how- 
ever, reduce this effect if the ensemble stars are roughly 
of equal colour and luminosity. We decided to correct for 
non-linearity for all sites except the Schmidt, which did not 
show measurable non-linear effects, and Laguna, where we 
only had measurements of the gain variations up to approx- 
imately 30000 ADU, which was signiflcantly lower than the 
intensity levels of most target stars. For LOAO, the data did 
not justify a description of the gain variation to be higher 
than a first-order polynomial fit, although a few points could 
indicate that higher-order features were present. The linear- 
ity calibration of the data from La Silla was based on a 
method described in Sect. lH.lTI hence no fit was applied to 
the data shown in Fig. |K| 

We obtained calibration images at La Silla for a new and 
more elegant method for determining the linearity proper- 
ties. This method provides a much more precise determina- 
tion of the CCD gain variations, which we will compare with 
the classical method in the next section. 



3.1.1 Ultra-high-precision method 

The basic concept of the method described in this section 
was first outlined by |3aldrv ( 19^ and [knudson (20o3) 
and w as developed into a fully applicable method bv lStelld 
i2002l) . Like the classical method, this method measures rel- 
ative variations in the CCD amplifier gain using flat-field 
images of different exposure times, but it differs by using 
flat fields that have a strong gradient e.g. by using a grism. 
Each fiat field showed a large smooth variation in light level 
from approximately the bias level to a significant fraction 
of the saturation limit, with the longest exposure reaching 
saturation (see Fig.|SJ. The advantage of this method is that 
the effect from instabilities in the light source used to ob- 
tain the flat fields is very small, because we are sampling 
a large range (in the longest exposures the entire range) of 
the CCD gain response in a single exposure. The resulting 
measurement precision of the gain variations is several or- 
ders of magnitude better than the classical method. Further, 
this method requires relatively few images to achieve high 
precision, making it very efflcient. 

We obtained spectral flats using the DFOSC spectro- 
graph on the Danish 1.54-m telescope (La Silla). Light vari- 
ation in one direction across the CCD was achieved using a 
grism to disperse the light from the slit illuminated by an 
internal telescope calibration lamp (Fig. (SJ. A series of im- 
ages were acquired in the following way; 3 x 30 s, 3 x 130 s, 
3 X 30 s, 3 X 250 s, 3 X 30 s, 3 X 370 s, 3 x 30 s, 3 x 390 s, 
3 X 30 s, 3 X 410 s, 3 x 30 s. The control exposures of 30 s 
enabled long-term drift in the fiat field-lamp to be removed. 
Although this improves the precision, it is not critical. Af- 
ter subtraction of overscan and bias, we corrected for the 
long-term drift of the fiat-field lamp and made an average 
(master) flat fleld for each exposure time. We then collapsed 
each master fiat field by averaging in one direction to obtain 
one dimensional intensity curves, as shown in Fig.Q 

Due to non-linear effects in the CCD, the intensity we 
measure in the nth row is: 

/(n) = rexp-S(n)-g(/(n)), (1) 
where Tcxp is the effective exposure time in seconds (after 
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Figure 7. Two intensity curves of different exposure times, 410 s 
and 130s (see text). The flat part of the 410s exposure at high 
row number is due to saturation of the CCD. 



correcting for dead-time of the shutter and short-term fluc- 
tuations in the hght source), g{I{n)) is the CCD amphfier 
gain as a function of the measured intensity and S defines 
the shape of the intensity curves, hence S{n) expresses the 
intensity in the nth row from a 1-s exposure. Although the 
actual gain variations are a function of the received flux, we 
assume the CCD amplifler has a well-defined output signal 
for every input signal. 

From two intensity curves with different exposure times, 
say 410 s and 130 s, we constructed the relative curve 



(2) 



/410b (n) _ 410s ff(74ios(n)) 
/i30s(n) 130 s g(/i30s(n))' 

where we have corrected for shutter dead-time and verified 
that the shape S was stable. In Fig. |S] (top panel) we show 
the curve /4i0s(n)//i30s(n) versus /410s (n), which we call a 
gain-ratio curve. At each intensity, Ji, this curve shows the 
gain ratio, g{Ii) / g{l2), where I2 = I\ ■ 130s/410s. Thus, 
finding g(I) for any intensity level requires inversion of the 
gain-ratio curve in an iterative process. 

We started out using a smoothed version of the mea- 
sured gain-ratio curve, which we denoted R\ (/), as a first es- 
timate for the actual underlying gain curve g{I) (see Fig.|Hl 
bottom panel). Then, assuming gi{I) = Ri{I), we computed 
a new gain-ratio curve, R2(I) ~ 5i(/)/ffi(/ ■ 130s/410s) 
for all /. The new estimate for the gain was then corrected 
by the relative deviation between Ri and R2 according to 
g2{I) ~ ^giil), etc. This iterative process stopped when 
Ri matched the measured Ri and the corresponding gi{I) 
was the desired gain curve. Before smoothing the measured 
gain-ratio curve (Fig. |5] top panel) we removed the noise at 
low intensities (0-1000 ADU) by replacing it with a linear 
fit to the data points from 1000 to 10000 ADU. In the exam- 
ple shown, the gain-ratio curve is very similar to the final 
gain curve. However, this is not a feature of the method, 
but is due to the gain characteristics of the particular CCD 
amplifier. 

To ensure that all features in the gain were detected, we 
examined gain-ratio curves based on different exposure-time 
ratios. If a feature, say a bump in the CCD gain, is periodic 
for increasing intensity and hence repeated at all pairs of 
intensities (7i,72) related as I2 — Ii ■ 410s/130s, it will 
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Figure 8. Top panel: Gain-ratio curve based on two inten- 
sity curves of different exposure time (410 s and 130 s). Bottom 
panel: Inversion from smoothed gain-ratio curve i?i (7) (dotted) 
to final gain-curve g{I) (solid). For comparison, the measurements 
from the classical method (Fig.|3 LaS) and their error bars (fx 
rms) are indicated. The uncertainty on g{I) is of the order 10 ppm. 
The smoothing of 7?i(7) introduces a systematic error in g(I) of 
roughly 1% near the saturation limit because the curve is steep. 



not show up in the gain-ratio curve based on 410 s and 130 s 
exposures (or any combination with the same exposure-time 
ratio). Our gain curves based on ffat fields with different 
exposure-time ratios, all showed an excellent match within 
the errors. 

In Fig. |H| (bottom panel) we compare our new method 
with results from the classical method for La Silla (Fig. (SJ. 
The two methods are in agreement with each other, but the 
series of flat-field images for our new method is significantly 
faster to obtain, provides the relative gain for all intensities 
and has a precision more than 100 times better. However, it 
requires temporally stable but spatially variable illumination 
of the CCD (e.g. spectral fiat fields) which is not possible at 
every telescope. 



4 ENSEMBLE PHOTOMETRY 

The goal of this project is to measure relative light varia- 
tions with very high precision. Hence, our approach is to 
obtain differential photometry taking advantage of the en- 
semble of stars in the FOV. Using an ensemble allows the 
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effects from atmospheric variations, common to all stars, to 
be removed from the time series. The number of stars in 
the ensemble ranged from 116 (in the small FOV of SAAO) 
to 3 58 (for LOAO). We used th e MOMF photometry pack- 
age iKieldsen fc Frandsen|[l993l to extract the photometry 
from the reduced images. It calculates differential photomet- 
ric time series by subtracting a reference time series which is 
a weighted average based on all ensemble stars. The weight 
given to each star is 1/rms"^, which ensures that stars with 
a high rms in their time series, such as faint stars and high- 
amplitude variables, are strongly suppressed. MOMF was de- 
veloped especially to produce time-series photometry from 
large numbers of images (in particular defocused images) of 
semi-crowded fields, similar to those obtained in this cam- 
paign. It combines PSF and aperture photometry. We chose 
10 stars, not necessarily the same for each site, to define 
the shape of the point-spread- function (PSF). These were 
all non-crowded bright stars, i.e. mostly red giant stars and 
a few bluer stars of the same luminosity (see Fig. . MOMF 
allows multiple apertures and calculates the total rms, crtotai, 
and the internal rms, aintemai, of the time series based on 
each aperture. The first is just the rms of the time series 
while the latter is calculated as 



2(iV 



N-l 



(3) 



where A'' is the number of points in the time series and rrii 
is the magnitude of the ith point in the series. For each site 
and each star we chose the aperture with the lowest crtotai 
in the time series. 




Figure 9. Cumulative distribution of d/a (solid line) for the en- 
tire data set of star No. 10. The dotted line is an analytical Gaus- 
sian distribution for comparison, and the dashed line indicates 
the threshold for our sigma clipping process. 
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5 IMPROVING THE PHOTOMETRY 
5.1 Iterative sigma clipping 

To improve the overall quality of the data we first removed 
outliers. We calculated the point-to-point deviation of each 
data point, i, relative to its neighbouring points as 



di^ rrii — 0.5 • (nii+i -I- rrn 



(4) 



where m is the magnitude. We then removed data points, 
rui, for which di > 3.0(7^, where o"; is the rms of d within 
a 3-hour interval around data point i. This was done in 
an iterative loop until no more points were removed, which 
converged after a few iterations. We illustrate the underlying 
statistics of this process in Fig. |5| which shows that the 
cumulative distribution of d/a has a significant non-white 
tail of deviating data points. Choosing the threshold is a 
trade-off between removing outliers and keeping statistically 
valid points. The threshold was chosen so that roughly 85% 
of the removed points were real outliers and 15% were valid 
points. We tested our sigma clipping method on generated 
random noise to verify that it was not too drastic in terms 
of removing extreme points from pure Gaussian noise. It 
removed less than 0.2%, corresponding to about 35 points 
of the time series from the entire campaign. This should 
be compared to the approximately 250 data points removed 
in total from the real time series for each star. Results of 
the sigma clipping on the real data were further verified by 
visual inspection of diagrams similar to what is shown in 
Fig.Hni 
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Figure 10. Photometric time series of star No. 10 for single nights 
from four sites. Outliers found by iterative sigma clipping are 
indicated with asterisks. The strong trend in the Kitt Peak data 
(curved time series) is due to colour-dependent extinction not 
removed by the ensemble normalisation of MOMF. The corrected 
time series, has been shifted upwards by 0.02 mag for clarity (see 
Sect.lCT. 



5.2 Colour extinction 

The data from Kitt Peak, which were the only ones ob- 
tained in the Johnson B filter, showed clear residual trends. 
This arises from an uncorrected wavelength dependence in 
the extinction (see Fig. 1101 . These trends were visible in 
the red giant stars because the reference time series calcu- 
lated by MOMF was dominated by bluer stars. Decorrelat- 
ing the differential photometry against airmass still left a 
lot of variation in the time series. Including more parame- 
ters, such as sky background, in the decorrelation process 
would affect the expected stellar signals significantly, which 
we verified with simulations. Subtracting a smoothed ver- 
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sion of the individual time series was also too harsh on the 
stellar signal. Instead we tried using only the red giant stars 
themselves in the ensemble. For some targets this provided 
good results but for most it did not. Adding the much fainter 
main-sequence stars of similar colour to increase the sample 
size did not improve the results. We finally chose to correct 
the colour term from the MOMF dif f erenti al photometry in 
a similar way to'Gillil and fc Brownl 1^8^. For each image 
we fitted a linear relation to the target stars: 



mfc 



(mfc) = ao+a-i- {B - V)k, 



(5) 



where is the magnitude of star k, [B — V)k is its colour 
and (rrifc) is the average of the time series. To correct star j, 
we subtracted a fit that did not include the star itself 



(m,)-(ao + ai-(B-V),). 



(6) 



This was to prevent stellar signal being removed by the pro- 
cess, which we have confirmed with simulations. There were 
12 red giant target stars with low noise levels, hence 11 stars 
were used in each fit, with a typical range of 1.00-1.25 mag 
m B — V colour. In Fig. 1101 (bottom right) we show the time 
series of one night of star No. 10 before and after correction 
of the colour term. The other sites showed weak effects from 
extinction, but these trends were not consistent from star 
to star or night by night, and the noise levels in the Fourier 
spectra did not improve if we performed the same correction 
as in the case of Kitt Peak. We therefore decided not to cor- 
rect for residual extinction in the differential photometry at 
any other site. 



(ii) We then adjusted the weights night by night to 
be consistent with the noise level (in amplitude), (Tamp, be- 
tween 300-900 /iHz in the Fou rier spe ctrum, requiring that 
o"amp E,1i = (Eq. 3 in iButler et al. 20041 The idea 
is that noise in this frequency range would have components 
that affect the noise at slightly lower frequency as well where 
we expect the stellar signal to be for the red giant stars. 
Choosing a frequency range within the expected range of 
the stellar signal would effectively down-weight stellar sig- 
nal, which is not desired. 

In Fig. 111! we plot our final estimates of Gi, including 
the adjustment multipliers for each night shown in the in- 
sets. The maximum adjustment was a factor of ~ 2. For 
some sites, the noise in the final Fourier spectra in the range 
300-900 /^Hz decreased by up to 20% after adjusting weights 
on a night-by-night basis, but in most cases it was a 5-10% 
decrease. We see that Gi vary significantly during the ob- 
serving run at many sites. For example, the range at Kitt 
Peak is 0.54-3.61 mmag (see Fig.im. 



6 ERROR BUDGET 

To establish whether the noise in the final time series was at 
the irreducible lower limit dominated by photon noise and 
atmospheric scintillation, we estimated each noise compo- 
nent and compared with the measured noise in the time 
series in a similar way as in prev ious investigations b y 
iGilliland fc Brownl (|l98i[l9923) and lGilhland et al.1 Jl993ll . 



5.3 Weight calculation 

To be able to detect solar-like oscillations in the red giant 
stars, it is crucial that we obtain noise levels as low as pos- 
sible in the frequency range where the oscillations are ex- 
pected to appear in the Fourier spectra of the time series. 
It is known that weighting time series of inhomogeneous 
data can signifi cantly improve the final signal-to-noise level 
/Handle il2003) . The important thing is that the final weights 
represent the true variance of the noise on time scales sim- 
ilar to the stellar signal one wants to d etect. We wil l use a 
weighting sch eme similar t o that used bv lButler et al.l (l2004l) 
and Kicldscn et alJ l|200^ to minimise the noise in ampli- 
tude, which includes the following two steps: 

(i) Calculate weights from the point-to-point variance 
[w, = l/a?). 

(ii) Adjust the weights to obtain agreement between the 
noise at the relevant frequencies in the Fourier spectrum and 
the weights as being represented by Wi — \/Gf. 

(i) The point-to-point variance was not supplied by 
the photometric reduction package and these values had 
to be estimated from the local variance of the time se- 
ries. We estimated the local scatter Gi (=\/ variance) for 
each data point i as the rms of the d array (Eq. |4j us- 
ing a moving boxcar. The width of the boxcar (5 hours) 
was chosen to minimise the noise (in amplitude) in the 
weighted Fourier spectrum. The spectrum was calculated as 
a weig hted discrete Fourier T ransform following the descrip- 
tion of lFrandsen et all (^9^. Having first removed outliers, 
we prevented good data from being down-weighted by bad 
neighbouring points during this process. 



Our total error budget comprised scintillation and 
counting statistics within the aperture; the latter included 
stellar light, sky background, level in the fiat field, and CCD 
read-out-noise. The contribution from scintillation was esti- 
mated as 



= 0.09D 



-2/3 3/2 



•X-'-exp(-V8000m)^,-;/^ (7) 

jKieldsen fc Frandsenlll99l E g. 3) using the factor of pro- 
portionalitv from T^un D the telescope diameter in 

centimetres, x the airmass, h the elevation of the telescope 
in metres and Texp the exposure time in seconds per im- 
age. For the counting statist ics we used the expression from 



age, l^or tne countmg statist ics we us 
iKieldsen fc FrandsenI lll99l Eq. 31) 



2 _ 21n2 1 

■^count TTTo ~^ 

W^TVes Csta 



2 Esky + CTpi^ON 
- ""''^AP 2 



(8) 



where W is the fuU-width-at-half-maximum of the stellar 
PSF in pixels, eg is the number of electrons per pixel in the 
fiat field, estar is the number of electrons from the star within 
the aperture, tap is the radius in pixels of the aperture, esky 
is the number of electrons per pixel in the sky background 
and (JRON is the CCD read-out-noise (per pixel) in electrons. 
Combining crscint and acount finally gives the estimated scat- 
ter 



_ / 2 2 N 

Cost — ^Tscint + (TcountJ 



1/2 



(9) 



The estimated scatter was dominated by scintillation for the 
brighter stars and by photon noise for the fainter stars. 
The magnitude at which the noise changed from being 
scintillation-dominated to photon-noise-dominated was in 
the range V = 10.5-12.0 mag but different from site to site. 
At a few sites the contribution to the counting statistics from 
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Figure 11. Scatter for each site for star No. 10. The insets show 
the multiphcation factor used to adjust ai for each night (see 
Sect. 1231 . The horizontal axes are the same as in the main panels. 



Figure 12. Measured scatter versus estimated scatter for the 
red giant stars. The measured scatter is (Tii^^crnal 

(Eq. EJ of the 

relative photometry based on the best night for each star. The 
estimated scatter is calculated from Eq. |2l which gives relative 
errors, and is scaled by a factor of 1.086 ppm//imag to put it on 
the magnitude scale. 
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the flat fleld was similar to the scintillation and hence sig- 
nificant for the brighter stars. In general, the read-out-noise 
and sky background could be neglected. 

In Table|21we give the measured internal scatter (Eq.|^ 
for each red giant star based on the full time series, which 
shows the overall quality of the data from star to star and 
from site to site. In general, star No. 10 had the lowest noise 
except for SAAO and LOAO. To compare with the estimated 
scatter (Eq.|5J we have, for each star and each site, measured 
the internal scatter (Eq.|3J for the best night, and the results 
are shown in Fig. 1121 There are other sources of noise not 
included in our error budget, which can explain why some 
stars fall significantly above the line of proportionality. Sat- 
uration of the CCD will increase the noise significantly. For 
several sites, stars No. 3, 4 and 11 were affected by satura- 
tion, which explains their higher noise levels. Close neigh- 
bouring stars can introduce higher noise in the photometry. 
Star No. 12 has three close neighbouring stars and at most 
sites it suffers from excess noise. If a star is located close to 
a bad column on the CCD the noise will also increase, which 
we see in some cases. 

In summary, for stars that are not affected by crowd- 
ing, saturation or bad columns, we generally see noise levels 
limited by photon and scintillation noise on the best nights. 
At one site (RCC) the noise is larger than estimated by a 
factor of two, which is unexplained. 

6.1 Noise comparison with previous campaigns 

The first campaign aimed at detecting solar-like oscilla- 
tions in M67 was carried out bv lGilliland fc Brownl Jl98d) . 
who used a 0.9-m telescope on two nights. Noise levels 
to ~ l.Smmag per minute integration we re attained for 
non-sa turated stars {mv > 12 mag). Later. iGilliland et alJ 
il99ll) observed M67 for two weeks from five sites using 0.6- 
m to 1.1-m class telescopes. The lowest rms in the time se- 
ries of the non-saturated stars (mv > 12.5) was 0.88 mmag 
per minute integration after high-pass filtering the data. For 
comparison, our best 0.6-m site (SOAO) showed an internal 
scatter (comparable to the rms of high-pass filtered data) 
down to ~ 1.6 mmag per minute integration, and the best 
1.0-m site (SAAO) showed ~ 1.2 mmag per minute inte- 
gration. A final, but direct, comparison can be made be- 
tween our Kitt Peak observatio ns of the red giant No. 10 
with those obtained bv Gillilan d fc Brownl (tl992a'l using the 
same telescope in a s imilar campaign. Star No. 10 (No. 7 in 
IGilliland etaDliggiTl was one of their targets with the low- 
est noise, which was 0.43 mmag per minute integration after 
high-pass filtering and decorrelating the data. In compari- 
son we obtained 0.50 mmag (based on the internal scatter, 
but without any high-pass filtering or decorrelation) . 

To summarise, we find noise levels as good as in previ- 
ous studies based on similar size telescopes. However, this 
campaign has provided significantly longer time series (six 
weeks compared to a maximum of two weeks) with better 
coverage than earlier comparable campaigns, which implies 
lower noise levels in the final Fourier spectra. 



6.2 Noise in the Fourier spectra 

The noise levels in the Fourier spectra obtained by our 
weighting scheme were 20-30 /imag for the red giant stars 
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Figure 13. Fourier spectrum of star No. 10 (in amplitude). The 
white line indicates the noise level in the range 300-900 ^Hz. The 
inset shows the spectral window, which is on the same frequency 
scale as the main panel. Each data point has been weighted ac- 
cording to the weighting scheme described in Sect. IF!^ 



that were not affected by excess noise from crowding, satu- 
ration or bad pixels. We therefore expect to be able to detect 
oscillations in the Fourier spectrum with S/N ^ 4. However, 
the detection threshold depends very mu ch on th e mode life- 
time, which is unknown for these stars JStello et al. 200^, 
and will require extensive simulations to quantify. This anal- 
ysis will be published in a subsequent paper. 

In Fig. 1131 we show the Fourier spectrum of star No. 10, 
which is one of the best, to illustrate the noise level we have 
obtained. The noise level, indicated with the white line, is 
the average in the range 300-900 /^Hz. The detailed pulsation 
analysis of all red giants will be presented by Stello et al. 
(in prep.). 



7 CONCLUSIONS 

We have collected 100 telescope nights of photometric multi- 
site data of the open stellar cluster M67 over a six-week 
period. The focus of this paper was the discussion of our ap- 
proach towards achieving the highest signal-to-noise ratio for 
the very low amplitude (50-500 /^mag) solar-like oscillations 
in the red giant stars. This included our careful reduction of 
the CCD images to obtain the lowest possible noise in the 
time-series data, and our weighting scheme to reduce the 
noise level in the Fourier spectrum. 

We have obtained a point-to-point scatter in the time- 
series photometry down to about 1 mmag for most sites, 
while the largest participating telescope reached 0.5 mmag 
(Fig. [TUl. These values are similar to those from pre- 
vious campaigns on M6 7 bv [C illiland fc Brown (198^, 
IGilliland etlD Jl99J) and lCilhland fc Brownl Il992bl) which 
all used telescopes of similar size but for shorter time spans 
(maximum two weeks) . Comparison of our best nights with 
known noise terms demonstrates that the attained point-to- 
point scatter is consistent with irreducible terms dominated 
by photon and scintillation noise for all sites but one (RCC), 
which shows extra noise of unknown origin (Fig. I12II . With 
these scatter values, our weighting scheme provided a mean 
noise level in the Fourier spectra of approximately 20 /imag 
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Table 2. Internal scatter tTinternal in mmag of the red giant stars (sorted by their luminosity). The internal 
scatter is based on the entire time series using Eq.|21 after ensemble normalisation, sigma clipping and, for 
Kitt Peak, correction for colour-dependent extinction (see Sect.|3. Star No. 10 has the lowest scatter at all 
sites except SAAO and LOAO. 



No. 


mv 


SSOi 


SSO2 


SOAO 


SAAO 


RCC 


Sch 


LaS 


LOAO 


Kitt 


Lag 


3 


9.72 


16.84 


23.69 


3.46 


- 


- 


13.11 


3.86 


8.18 


- 


7.75 


11 


9.69 


26.76 


10.28 


3.38 


4.11 


- 


16.81 


1.73 


10.10 


- 


2.88 


4 


10.30 


2.37 


9.55 


1.87 


1.84 


5.16 


2.83 


1.39 


3.33 


2.77 


2.94 


21 


10.47 


- 


- 


20.63 


- 


- 


- 


- 


3.77 


- 


- 


8 


10.48 


2.37 


8.13 


2.00 


1.70 


3.96 


2.58 


1.38 


2.95 


1.46 


2.51 


9 


10.48 


2.73 


9.44 


1.81 






3.16 


1.42 


2.94 


1.22 


3.24 


10 


10.55 


2.35 


2.02 


1.44 


2.04 


3.87 


2.31 


1.36 


3.09 


1.11 


1.95 


20 


10.55 


2.86 


17.43 


2.85 






2.93 


1.58 


3.13 




5.57 


2 


10.59 


2.54 


4.83 


2.09 






2.91 


1.63 


2.92 


1.28 


4.09 


18 


10.58 


3.24 


16.04 


4.62 






3.12 


1.41 


3.07 




7.64 


16 


10.76 


2.96 


9.69 


3.16 






3.28 


1.55 


7.93 




5.15 


5 


11.20 


2.44 


5.22 


3.04 


1.82 


4.68 


3.68 


1.96 


2.63 


1.62 


2.78 


17 


11.33 


2.72 


5.25 


2.47 




5.66 


3.22 


3.50 


2.45 


1.91 


3.20 


7 


11.44 


2.83 


13.65 


2.83 


2.00 


6.73 


3.31 


1.54 


2.57 


1.61 


17.38 


19 


11.52 


2.91 


4.86 


2.85 






3.20 


1.64 


2.95 


4.81 


3.60 


15 


11.63 


2.50 


4.54 


3.42 




5.33 


3.54 


6.22 


2.66 


1.83 


2.88 


14 


12.09 


3.02 


5.40 


3.96 




4.75 


4.22 


1.66 


3.23 


3.59 


2.90 


12 


12.11 


5.39 


8.88 


8.47 


3.98 


8.19 


4.30 


8.36 


6.04 




9.02 


13 


12.23 


3.49 


8.74 


5.41 






4.75 


3.72 


3.32 


3.53 


3.26 


6 


12.31 




13.16 


5.02 






6.24 




3.47 







(in amplitude), which would allow us to detect solar-like 
oscillations in the red giant stars with S/N ^ 4 assuming 
L/M-scaling. 
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